Low work function of the (1000) Ca 2 N surface 
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Polymer diodes require cathodes that do not corrode the polymer but do have low work function 
to minimize the electron injection barrier. First-principles calculations demonstrate that the work 
function of the (1000) surface of the compound Ca2N is half an eV lower than that of the elemental 
metal Ca (2.35 vs. 2.87 eV). Moreover its reactivity is expected to be smaller. This makes Ca2N 
an interesting candidate to replace calcium as cathode material for polymer light emitting diode 
devices. [PACS: 73.30 +y] 
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One of the great challenges for polymer light emitting 
diodes is the electron injection barrier as well as the per- 
formance degradation caused by chemical reactions of 
the cathode with the polymer. 1 State-of-the-art devices 2 
use PPV (poly phenylene vinylene) as electro lumines- 
cence material. Often the cathode is made out of cal- 
cium, because of its low work function (2.87 eV) 3 and 
the presumed alignment of its Fermi level with the low- 
est unoccupied molecular orbital of PPV (electron affin- 
ity 2.73-2.8 eV). 4,5 Biases actually used are much higher 
than the optical band gap of PPV (2.4 eV) 6 namely in 
the order of 10 V, 1 which suggests that there is some 
barrier formation between the cathode and the poly- 
mer. This is further confirmed by the observation that 
the decrease of performance is much less for devices in 
which the cathode is formed at a small residual oxygen 
pressure. 7 In such a device some oxidation prevents re- 
action of the cathode with the polymer. 

It would be a major breakthrough when stable, low 
work function metals are found. Sub-nitrides are fas- 
cinating candidates because, as in the case of cesium 
sub-oxides, 8 ' 9 it is expected that through quantum con- 
finement their work functions are lower than those of the 
elemental metals while at the same time ionic bonding 
reduces reactivity. 

The Ca2N crystal 10 and electronic structure (bulk and 
single slab) 11 are known. It is an ionic compound, built 
out of slabs of alternating (Ca-N-Ca) hexagonal layers. 
The interslab distance between Ca layers is 3.81 A and 
Ca layers inside a slab are only 2.45 A apart. This is 
much smaller than the ordinary fee Ca (111) layer dis- 
tance (3.18 A) and can be attributed to the ionic binding. 
The band structure shows fully occupied N 2s and N 2p 
states and a quasi 2D free-electron state in the space be- 
tween the slabs which further confirms the idea that this 
is an ionic compound. 

In this paper the Ca2N work function ($) and surface 
energy (E$) are determined for both the (0001) and the 
(1000) surface. An indication of the materials stability 
is obtained from its binding energy and degree of surface 
relaxation. 
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Fig. 1. 

aged over (0001) planes and potentials (relative to Ef, eV) 
averaged over (0001) planes and bulk unit cells (solid/dotted 
lines respectively) as function of position (A). The horizontal 
line at -3.83 eV refers to the average bulk potential. 

First-principles calculations were carried out based 
on density functional theory in the local den- 
sity approximation 12,13 with generalized gradient 
corrections. 14 The total energy and molecular dy- 
namics program VASP (Vienna Ab-initio Simulation 
Package) 15,16 was used which has the projector aug- 
mented wave method 17,18 implemented. Non-linear core 
corrections 19 were applied for calcium. The unit cells 
contained 18 and 42 atoms for the (0001) and the (1000) 
surfaces respectively. The Kohn-Sham orbitals were ex- 
panded in plane waves with cutoffs of 37 Ry. 12 • 12 ■ 1 
and 1-12-4 Monkhorst-Pack 21 k-point grids respectively 
were used to sample the Brillouin zones resulting in 74 
and 26 k-points respectively in their irreducible parts. 

$ is defined as the minimum amount of energy it costs 
to extract an electron from a metal. The fact that it de- 
pends on the type of surface, is somewhat counter intu- 
itive, it seems to contradict the conservation of energy. 
However, the differences in work function are compen- 
sated by differences in kinetic energy of the electron at 
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large distances. 
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Fig. 2. Atomic positions after relaxation of the (1000) 
surface as seen from the (1206) direction. Dark circles are 
nitrogens and light ones calcium. 

In order to determine $, two numbers are needed: 
the maximum energy of an electron inside the material, 
defined to be the Fermi level (Ep), and the minimum 
energy outside the material, the vacuum potential. Be- 
cause in a finite slab calculation Ep cannot be accurately 
determined, this number is first obtained from a bulk 
calculation. Then for the finite slab the average poten- 
tial of the middle layers (which is an accurate quantity) 
must be aligned with that of the bulk. 23 In this way the 
work function for the (0001) surface is determined us- 
ing 6 slabs of Ca-N-Ca and 19 A of vacuum. 24 It can be 
extracted from Fig. 1 and equals 3.43 eV. 

Another unit cell was made to find <E> for the (1000) 
surface of Ca2N. It contained 14 layers and 20 A of vac- 
uum. At this surface several bonds are severed and it 
is necessary to relax the atomic positions. The result is 
shown in Fig. 2 where the view is along the displacement 
vector from one slab to the next. The nitrogen atoms at 
the surface are displaced a little (0.24 A) inwards, the 
nitrogens in the second layer a little (0.17 A) outwards, 
while the surface calciums tend to move in the direction 
of the surface nitrogens. This can be understood from an 
ionic point of view. Nitrogen favors an environment of 
high electron density while calcium shows the opposite 
trend. Although the structure has changed very little, 
the effect of the relaxation on $ is substantial (Fig. 3). 
Before relaxation we find a work function of 2.56 eV, 
already lower than that of pure calcium, but after re- 
laxation it has even decreased to 2.35 eV. Not only is 
this comparable to the smallest work function for an el- 
ement (Cs, 2.14 eV), 3 but also the difference of 1.08 eV 
with the (0001) surface is larger than that between any 



two surfaces of tungsten, 25 which is known for its large 
work function anisotropy. Both the drop in <J> after relax- 
ation and the large surfacc-anisotropy can be explained 
by the Smoluchowski model 26 and not by quantum con- 
finement as this is a bulk model that does not say any- 
thing about anisotropies and/or relaxation. According 
to Smoluchowski on the other hand, the more "open" a 
surface is the more the charge density can be smoothed, 
and this lowers However more-open surfaces tend to 
have higher surface energies 27 and this makes them less 
stable. For the (0001) and (relaxed) (1000) surface we 
find Eg is 2.7 and 5.3 eV/nm 2 respectively which con- 
firms the trend. 




Fig. 3. Charge densities averaged over (1000) planes 
(arb. units) and potentials (relative to Ef, eV) averaged 
over (1000) planes and bulk unit cells, for both the relaxed 
(solid lines) and the unrelaxed (dashed lines) structure as 
function of position (A). The potentials converge to 2.35 eV 
and 2.56 eV in the vacuum respectively. The horizontal line 
is the average bulk potential. 



Because in Rcf. 10 it is suggested Ca2N can also be 
made with a larger c axis (interslab separation), we 
examined whether this changes the (1000) work func- 
tion. After enlarging c with 2% and subsequent re- 
laxation of the atoms, the structural changes are small 
and the trends already indicated are a little more pro- 
nounced. The slab thickness however has decreased to 
almost (+0.5%) its former value. The averaged charge 
density and local potential are plotted in Fig. 4. Be- 
fore relaxation <!> is decreased from 2.35 to 2.24 eV, but 
relaxation increases it again to 2.35 eV. Again this is 
in accordance with the Smoluchowski trend because a 
thicker slab makes the surface more open. We conclude 
that the work function does not (significantly) change 
when c is enlarged. 
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Fig. 4. Charge densities averaged over (1000) planes 
(arb. units) and potentials (relative to Ef, eV) averaged 
over (1000) planes and bulk unit cells, for both the relaxed 
(solid lines) and the unrelaxed (dashed lines) structure for a 
2% larger c axis as function of position (A). The potentials 
converge to 2.35 eV and 2.24 eV in the vacuum respectively. 
The horizontal line is the average bulk potential and it can 
be seen that relaxation has no observable effect on the charge 
density. 

Finally the reactivity of Ca2N is briefly discussed. In- 
spection of the binding energy shows that it is 5.6 eV per 
formula unit (with respect to half a nitrogen molecule 
and 2 calcium atoms) for a single slab and 0.6 eV more 
for the bulk. This gives evidence that the reactivity of 
Ca2N is smaller than that of calcium, of which the bind- 
ing energy is just 1.9 eV per atom. It was shown earlier 
that relaxation of the atomic positions after bonds are 
cut or c is enlarged, does not have a substantial effect on 
the structure. These facts taken together, give a indi- 
cation that the reactivity of this sub-nitride is less than 
that of pure calcium, although it is still air sensitive. 10 

In summary, using first- principles calculations, we 
have shown that the work function of Ca 2 N has an 
anisotropy of more than 1 eV, with a minimum of 2.35 eV 
(cf. Ca, $ = 2.87 eV) for the (1000) surface. Its sur- 
face energy (5.3 eV/nm 2 ) however is almost the double 
of that of the (0001) surface. The low work function 
does not (significantly) depend on the interslab separa- 
tion and we argued that this material is not as reactive 
as calcium. The lower work function and reduced reac- 
tivity make Ca 2 N a promising metal to replace Ca in 
cathodes of polymer light emitting diodes. 
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